home *** CD-ROM | disk | FTP | other *** search
/ CU Amiga Super CD-ROM 18 / CU Amiga Magazine's Super CD-ROM 18 (1997)(EMAP Images)(GB)[!][issue 1998-01].iso / CUCD / Utilities / Laplace / Projects / physic_PI_41 < prev    next >
Text File  |  1996-06-28  |  691b  |  45 lines
  1. #LAPLACE
  2. #USEFLOAT=FALSE
  3. #SIMPLIFY=TRUE
  4. #ERRBND=15
  5. #EXPDSP=5
  6. #PREC=2
  7. #F1=f(x)=
  8. #F2=g(x)=
  9. #F3=h(x)=f(x)/g(x)
  10. #F4=matrix(
  11. #F5=vector(
  12. #F6=diff(f)
  13. #F7=diff(g)
  14. #F8=diff(h)
  15. #F9=taylor(f,0)
  16. #F10=taylor(f,1)
  17. #F11=%(run type s:startup-sequence)
  18. #F12=
  19. #F13=
  20. #F14=
  21. #F15=
  22. #F16=
  23. #F17=
  24. #F18=
  25. #F19=
  26. #F20=
  27. H=set(72.8-124,111-124,129.2-124,155-124,180.5-124,204.5-124,229.8-124,254.5-124,278-124,300.5-124,328-124)
  28. h_0=-256/5
  29. alpha=0.003674
  30. b=775.87
  31. t_z=19.38
  32. t(h)=(h-h_0)/(alpha*(h_0+b))
  33. f(t)=1+0.008*(t-t_z)/(t_z+273)
  34. T(h)=((h+b)*f(t(h)) - (h_0+b)*f(0))/(alpha*(h_0+b)*f(0))
  35. t(204)
  36. T(-51.2)
  37.  
  38. c_w=4.182,m_w=205.3
  39. t_1=21.5,t_2=51
  40. W(t)=c_w*m_w*(t_2-t)/(t-t_1)
  41. dt=0.2
  42. W'(t)=eval(diff(W,t))
  43. abs(W'(49.1))*0.2
  44.  
  45.